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Abstract. This work's purpose is to understand the dynamics of some social systems 
whose properties can be captured by certain iterated function systems. To achieve this 
intension, we start from the theory of iterated function systems, and then we study two 
specific economic models on random utility function and optimal stochastic growth. 



1. Iterated Function Systems 

1.1. Discrete Dynamical Systems. The transition rule in a discrete dynamical system 
can be determined solely by a single function mapping the state space into itself. Let the 
state space X be a metric space with the metric d : X x X ^ M+ , and let the time domain 
be T = Z. 

Definition 1.1. A discrete dynamical system on X is a pair {X, /), where f : X ^ X 
is a continuous function, such that Xn+i = f{xn) for all Xn-,Xn+i S X and all n G Z. 

If we consider a state x, by the iterated transition rule, we can have a forward series 
of states as f{x), f{f{x)), f{f{f{x))), and so on, and also a backward series of states 
as f~^{f~^{x)), f~^{f~^if~^ix))), and so on. To make our notations concise 

enough, we define in a recursive way, say = /""^ o / for n G Z, and = idx, where 
idx is the identity function on X. We have a discrete motion passing through x being 

...,/-"(x),...,rHx),x,/(x),...,r(x),..., 

where the transition rule between the time m with a state x £ X and the time m + n is 
/"(x), for any m,n £7^. 

The trajectory passing through x is 

and similarly the positive and negative semi-trajectories are 

7"'(x)= U {/"(^)}' and7-(x)= [j {/"(x)}. 

We know that 7^(x) U 7~(x) = 7(x), and 7'''(x) n 7~(x) 3 {x}. If 7(x) = {x} or 
equivalently /(x) = x, we call x is an equilibrium state, or in terms of functions, a fixed 
state. 

A set of states S X is said to be invariant if f{S) = S, where f{S) := {/(x) : x G S}. 
The cj-limit set of a state x is w(x), which contains all the cj-limit states of x. If uj{x) ^ 0, 
then it must be invariant, say /(cij(x)) = uj{x). 
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A state X is said to be periodic, if there exists a A; G such that f^{x) = x. The 
period p of 7'^(x) for any periodic state x is defined as the minimal integer k such that 
f{x) = X. 

If the period of 7^(a;) is 1, then f{x) = x, so x is a fixed state. If the period of 7+(x) 
is p < +00, then 'y'^{x) = {x, f{x), . . . , f^{x)}. A state x is said to be finally periodic, if 
there exists an m G Z+ such that /"(x) is a periodic state for all n > m, or equivalently 
/"+P(a;) = /"(x) for all n > m, where p is the period, x is said to be asymptotically 
periodic, if there exists a y e X such that lim„_^+oo d{f'^{x), f"'{y)) = 0. 

A set of states A C X is called an attractor, if there exists a neighborhood N{A, e) of 
A such that f{N{A,e)) C N{A,e), and 

u;{NiA,e))= f] r{NiA,e)) = A. 

Now consider a compact state space AT C M, we can state the following theorem that 
any periodic patterns can be implied by the existence of a period 3 trajectory. 

Theorem 1.1. If / : X — t- X is continuous where X is an interval in M, and there is 
a trajectory of period 3, namely there are three states x,y,z G X such that /(x) = y, 
f{y) = z, and f{z) = x, then there exist trajectories of period n, for any n G N. 

Proof. See the proof of Theorem 10.1 in Devaney [De]. □ 

In general, this result is a special case of Sarkovskii's theorem. Define the Sarkovskii's 
ordering of all the natural numbers by >s, if there is a trajectory of period n, then there 
also exist trajectories of period m, where n >$ m. Since 3 is the first natural number in 
Sarkovskii's ordering, and all the powers of 2 in the order 

■ • • >5 2" >5 2"-l >s • • ■ >5 22 >5 2 >5 1 

are the last group of numbers in Sarkovskii's ordering, a period 3 trajectory implies all 
the other possible periods, and if we have a period not a power of 2, then we should have 
all the periods of the powers of 2. 

The discrete dynamical system (A, /) is chaotic, if 

(i) for any x G A and e > 0, there exist 5 > and n G Z+ such that rf(/"'(x), > 

e for all y G A(x, 5), 

(ii) for any pair of open sets of states and 52, there exists an n G such that 

r (Si) n S2 7^ 0. 

So a chaotic discrete dynamical system means that its dynamics is sensitively dependent 
on the initial states, and its states are transitive. 

We also want to give a definition of chaos in the sense of Li and Yorke [LY], in which 
A C M is compact, and / : A — >■ A is of class C^. 

Definition 1.2. A discrete dynamical system (A,/) is said to be non-periodically 
chaotic, if there exists an uncountable set S C A such that 

(i) limsup„^+ooC^(/"(a;),/"(y)) > 0, and liminf„_^+oo /"(y)) = for all 
x,y € S with x ^ y, 

(ii) for all x,z & S with z is periodic, limsup„_j._|_^ d(/"'(x), /"(z)) > 0. 



ITERATED FUNCTION SYSTEMS WITH ECONOMIC APPLICATIONS 



3 



This definition of chaos is sHghtly weaker than the definition we give before. If a 

discrete dynamical system on X C ]R is chaotic, then it must be also non-periodically 
chaotic. On the other hand, if a discrete dynamical system is non-periodically chaotic, 
it is not necessary to be chaotic. Notice that the condition (ii) in the definition of non- 
periodic chaos means that there is no asymptotically periodic state in X. Suppose u is a 
state not asymptotically periodic, so u;{u) has infinitely many states, and there exists a 
unique minimal invariant set V C Ld{u), such that V = io(v) having infinitely many states, 
where v e X. Consider the set of states U = X\V, then /"(F) n [/" = for all n G Z+, 
so it is not transitive, which hence means that it is not chaotic. 

1.2. Iterated Function Systems. In this section, we consider a set of contractive func- 
tions defined on the state space X with a metric d : X x X ^ IR+, in which a function 
f : X ^ X is said to be contractive if there exists a A G (0, 1) such that d(^f{x), f{y)) < 
Xd{x,y) for all x,y e X. In the following paragraphs, we denote 7jv as an index set with 
N elements, where AT > 2 is finite. 

Definition 1.3. Let F be a collection of > 2 continuous contractive functions, say 
F = {fi : i £ In} where fi : X ^ X. Then the pair {X,F) is said to be an iterated 
function system, if {X,fi) is a discrete dynamical system for each i G In- 

Suppose {X, d) is compact, and Q{X) is the set of all nonempty compact subsets of X. 
For any U,V ^ Q{^)i the Hausdorff metric between U and V is 

dH{U,V)= sup {d{u,V),d{v,U)}, 

ueu,vev 

where d{x, Y) := inf^gy d{x, y) is the metric between a state x £ X and a nonempty set 
Y C X. Then Q{X) with the Hausdorff metric dn ■ QiX) x Q{X) M+ is a compact 
metric space. 

Consider a mapping H : Q{X) — >■ Q{X) such that for any B G Q{X), we have 

H{B) = (j MB), 

where H is called a Hutchinson operator. Define i?" in a recursive way, say H"- = H"^~^oH 
for 77. G Z, and = idQ(x)) where again idgf^x) is an identity mapping on Q(X). 

Definition 1.4. A G Q{X) is called an attractor of the iterated function system {X, F), 
if H{A) = A, and there exists a neighborhood N{A,e) G Q{X) such that H[N{A,£)) C 
N{A,e) and ^£2+ H-{N{A,e)) = A. 

Obviously, if A is an attractor of {X,F), then there exists a neighborhood N{A,e) G 
Q{X) such that lim„^+oo ^"(^(^, e)) = A. The largest such neighborhood is called the 
basin of the attractor A, denoted as B{A). 

Proposition 1.1. The mapping H : Q{X) — Q(X) has a unique fixed point A G Q{X), 
i.e. H{A) = A, and hence the iterated function system {X,F) has a unique attractor A. 

Proof. Since /j is a contractive function, there exists a A^ such that d(^fi{x),fi{y)) < 
Xid{x,y) for all x,y E X, and any i G In- Define A = maxjg/^ Aj, then 

dH{H{U),H{V)) < sup dH{fi{U),f,{V)) < sup {X.,dH{U,V)} < Xdn{U,V), 
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for all U,V E Q{X). So H has a unique fixed point, say A, which is also the attractor of 
iX,F). □ 

Now consider the space := In x In x • • • , and for any // G I^, we write = 
(^i)A*2)''' il^ki''')^ where fi^ € In for all A; G N. For any £ I^, define the Baire 
metric as dsi/J', v) = where j is the first index when / Vk- Now we have a compact 
metric space (7^,^^). For any given set of states S C X, we can define a mapping 
C -.1^^ Q{X) such that 

keN 

where n defines a trajectory starting from any state x € as 

7+(x) = {UAx)} U {f,, o f^^(x)} u • • • U o o . . . o f^^{x)} U • • • . 

If 5 C ^(A), then for any jj, G we have C{ii,S) C ^4, and especially we have 
C{I'^ ,B{A)) = A. So we claim that the attractor of an iterated function system can be 
achieved by all the possible trajectories determined by /i G Ij^, in which the transition 
function is f^^, at time /c, where fi^ G In- 

We can also impose a probability measure on such that ji can be realized with some 
probability density. In a simplified situation, the probability measure over the whole space 
could be represented by a stationary discrete probability measure tt over In such that 
^ie/jv ^^^-^ ~ ^^^^^ independent with the time and the state. At any time G Z, 
the transition rule is chosen at random from F according to the probability measure tt 
over In- 

Definition 1.5. Let {X,F) be an iterated function system, where F = {/j : z G In} 

where fi : X ^ X. If tt is probability measure on In such that X^ie/^v ^^^^ ~ ^^^^ 
triplet (X, F, tt) is said to be an iterated random function system. 

We use the random variable to represent the index chosen at time k, and the related 
random transition rule is fa,, at time k, so we have the random iteration 

Zk+l = fak+ii^k), 

where Zj- is the random state with the support of X for all k £N. Being similar with the 
above notations, a is defined as (cri, (72, . . . , CTfc, . . .). 

If the initial state is x, its random trajectory is written as T~^{x) = UfceNi'^fe}' where 
Zi = faAx), and = fa^.{Zk-i) for G N \ {!}. Clearly, the stochastic process {Zk)kGfi 
is a Markov chain. Given a state a:; G X at time k, the next state located in y C X has a 
probability 

P(x,y)= ^7r(z)ly(/i(x)), 

where 

1, if yeY 
0, ify^Y- 

For any Borel subsets Y of X, there exists an invariant probability measure p such that 



iY(y) 

, there ' 

p(Y)= [ P{x,Y)dp{x)=Y,'^ii)pUr\Y))- 



ieiN 
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Such a probability measure is called the 7r-balanced measure on {X,F,tt), as discussed in 
Barnsley and Dcmko [BD]. The support of p, defined as R{p) = {x £ X : p[x) ^ 0}, 
satisfies R{p) = [Ji^ij^ fi(R{p))^ a^^d hence R{p) = Ahy Proposition 1.1, which means 
that the support of a 7r-balanced measure on (X, F, vr) is the unique attractor A of (X, F) 
for any possible tt. 

Therefore, we claim that the attractor in (X, F) can be achieved by a random itera- 
tion generated from [X, F, tt) after a sufficiently long time. For any given state x G X, 
Pl^gjsj r"'"(Zfc) = A with probability 1, where G r+(x) for all fc G N, which means that 
there exists a random trajectory starting from some states in another random trajectory 
whose initial state is given at the beginning, and this random trajectory will cover the 
attractor eventually. Or in terms of limit states, the random oj-limit set of any given initial 
state X, say will be equal to A for sure. 

Now consider a special case in which the state space X is a subset of the n-dimensional 
Euclidean space E", and /j : X — >• X is affine, say fi{x) = DiX + Cj, where x G X is an 
n X 1 vector, in an n x n matrix, and Cj is an n x 1 vector, where i G In- Suppose again 
there is a probability measure vr over /jv, then the random iteration is Z^+i = DZ^ + e, 
where D and e depends on the randomly chosen index G In according to vr. We know 
that (^fc)fcGN is a Markov chain, so this random iteration can be equivalently represented 
as a general stochastic process Zj^^i = D^+iZ^ + 6^+1, where Z^ is the random state with 
the support of X, and D}^,ek are i.i.d. random matrices. 

Example 1.1. Suppose X = [0,1], fa{x) = x/3, h{x) = ,r/3 + 2/3, F = {fajb], 
h = {o-jb}, and 7r(a) = 7r(6) = 0.5. This iterated random function system is equivalent 
with the following autoregressive process 

Zk+i = Zk/3 + Ek+i, k G Z+, 

where Zq = x G [0, 1] is the given initial state, and is an i.i.d. random variable taking 
and 2/3 with equal probability. (X, F) has an attractor 

oo 

^={Y1 ^^/^* • ^ -fO' Vi G n}, 

i=l 

which is actually the Cantor ternary set. 

Let Bn := {EZn^i/^' ■ Xi G {0,2}, G N}, then A = Bi, and faiA) = A/3 = 
B2 = B2 + 0/3, h{A) = A/3 + 2/3 = B2 + 2/3, where A/3 = {a/3 : a e A}, and 
B2 + c:={b + c:b£B2}forb = 0, 2/3, so fa{A) U fb{A) = Bi = A, which confirms that 
A is the unique attractor of (X, F) . 

The attractor A is the support of a vr-balanced measure on (X, tt), and thus it can 
be achieved by a random trajectory r+(.7;) starting from any x G [0,1]. The random 
trajectory can be represented by the associated autoregressive process. 

If the initial state x e A, then we have r+(x) = A for sure. If x G [0, 1] \ A, then 
there exists xi, X2, ■ ■ ■ , Xm S {0, 2} such that x = ^ Xj/3* + m{x), where m{x) < 1/3™". 
If m is sufficiently large, then m{x) is close enough to 0, while if m — t- +00, m{x) — )■ 0. 
Zi = x+ei = y+m{x)/3, where y = ^™iXi/3*+ei G A, and generally, Zk = z+m{x)/3^ , 
where again z G A, so after a finite period of time, say Z, we will have Zi G A, and hence 
we have T^{Zi) = A with probability 1. Thus for any x G [0, 1], this stochastic process 
will cover A at last, say ^l{x) = f^i^^fqT^{Zk) = A, where Z^ G T~^[x). □ 
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1.3. Strange Attractors. The concept of strange attractors first appeared in a paper 

on turbulence by Ruelle and Takens [RT], but was not defined precisely. A definition of 
strange attractors later was given in a popular article by Ruelle [Ru], but according to 
that definition, the attractors of any chaotic dynamical systems are strange. We want to 
restrict this definition by stating that an attractor is strange if the dynamical system is 
chaotic, and it should a fractal. The fractals are defined by Mandelbrot as 

Definition 1.6. A set is said to be a fractal, if its Hausdorff-Besicovitch dimension 

strictly exceeds its topological dimension. 

Let {X, d) be a metric space as usual, and Z C X, the topological dimension of Z 
denoted as dimi(Z) is an integer belongs to Z_|_LJ{— 1, +oo}. If Z = 0, then dimj;,(Z) = —1. 
If Z is a discrete set, then dim/,(Z) = 0. In general, dmiL^Z) = n, if for any open cover 
of Z, any state z £ Z can be contained in n + 1 open sets of that cover at most. If there 
does not exist such a finite n, then dimi(Z) = +oo. 

For any d > 0, we can define the ^-approximate d-dimensional Hausdorff measure of Z 

as 

nl{Z) = inf { {di&m{Ui)y -.[jUiDZ, diam(C/i) < 

i£l iel 

where diam(C/j) := sup{d{zi, Z2) ■ zi,Z2 G Ui}, / is a countable index set, and {Ui}iQi is 
an open 5-cover of Z. The d-dimensional Hausdorff measure of Z is defined as n'^{Z) = 
lim5_=.o HgiZ). Then the Hausdorff-Besicovitch dimension of Z is defined as dimniZ) such 
that 

dimH(Z) = mt{d e R+ : W^iZ) = 0} = sup{d G M+ : W^iZ) = 00}. 

In practice, when we need to determine the dimension of a fractal, the more useful 
notion is the box-counting dimension. The upper and lower box-counting dimensions of 

Z C. X are defined as 

logiV(Z,e) ,^,_,,^,AogN{Z,e) 



dimB(Z) = limsup — j — ' , and dim p(Z) = liminf 



,^0" log(l/e) ' ^ .-.0 log(l/£) ' 

where N{Z, e) in the minimal number of sets as required for an open £-cover of Z. As we 
suppose N{Z, e) ~ A(l/e)'^ by extending the observations of regular cases in the Euclidean 
space, we have \ogN{Z, e) = log A-|-dlog(l/£), where A is a positive constant. That's why 
we define the box-counting dimension by the limit of the ration between log N{Z, e) and 
log(l/£). 

Once we have dimB(Z) = dim^(Z) for a set of states Z, we say there exists the box- 
counting dimension, which is denoted as 



d\m.B{Z) = dmiB{Z) = dim ^fZ). 



In general, we have dmiH{Z) < dim p(Z) < dims(Z), and dim//(Z) < dims(Z) if 
dim p(Z) = d\m.B{Z). For the unique attractor A of an iterated function system {X,F), 
if X C E" (or in some other special settings), we have d\mH{A) = d\mB{A), and hence 
we can calculate the box-counting dimension of A to get exactly its Hausdorff-Besicovitch 
dimension. If A is strange, we have dim.H{A) = dims(^) > d\m.L{A). 

Consider an iterated function system (X, F) has a set of affine functions with similarities 
/j : R" ^ M" for all ieiN- We know Uig/^^ •^»(^) = ^' ^^^^^ ^ attractor of {X, F). 

If we assume /i(^)'s are not overlapping, which means that there exists an open subset 
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D oi A such that fi{D) n fj{D) = for all i,j e In with i ^ j, and \J-^j^ fi{D) C D, 
then the similarity dimension of A, which is equal to dimB(^) and dimH{A), is the unique 
solution of the equation 

ieiN 

where Aj < 1 is the Lipschitz constant of the function /j for all i £ In- 

The intuitive proof of the above equation is clear. Suppose N{A, e) is again the minimal 
number of sets required to fully e-cover A, then N{A,e) X{l/sf- Notice that A = 
UieiN fi^^")' and /i(^)'s are not overlapping, so we have N{A,e) = Y^^^j^ N{fi{A),e). 
But iV(/i(^),e) = N{A,£/Xi) as A is scaled by the Lipschitz constant Aj of fi{A). Then 
N{A,e) = Eie7jv^(^'^/^i)' '^^d hence (1/e)'^ = Eie/^C^i/e)'^ when e ^ 0, so we get 

In the following paragraphs, we will consider some examples. At first we want to revisit 
the Cantor set that has been discussed in Example 1.1. 

Definition 1.7. Let Di = {0, 1} for all i G N with the discrete topology, and write 
{0, 1}'^ = HieN-^i' ^^^^ any set homomorphic to {0, 1}'^ is called a Cantor set. 

Example 1.2. Let X = [0,1], and define a real number c G (2, +oo). We remove the 
interval (1/c, 1 — 1/c) from [0, 1] having 

Ci = [0, 1/c] U [1 - 1/c, 1], 

and we delete the (1 — 2/c) ratio of middle parts in [0, 1/c] and [1 — 1/c, 1] respectively 
obtaining 

C2 = [0, l/c^] U [1/c - l/c^l/c] U [1 - 1/c, 1 - 1/c+l/c^] U [1 - l/c2,l]. 

In general, we have Ck+i = Ck/c U (Ck/c + (1 — 1/c)). When k — )• +00, we have the 
Cantor set 

C= [^Xi/d : Xi e {0,c- 1}, Vi G n}, 

which is the attractor of (X, F) with /a(x) = x/c and fb{x) = x/c + (1 — 1/c). 

If taking e = l/c^, we have N{C,e) = 2'^ as we need 2*^ open sets to e-cover C 
with 2^ disconnected intervals. By the definition of the box-counting dimension, we have 
dimij(C) = dimB(C) = limjt^oo log 2^/ log c'' = log 2/ log c. 

Alternatively, since C is the attractor of an affine function systems with similarities, 
and we have their scaling numbers Aa = Aj, = 1/c, so dimB(C) is the unique solution of the 
equation (l/c)"* + (1/c)'^ = 1. Thus dimB(C) = logc2 = log 2/ log c. Since c G (2, +00), 
dims(C) G (0,1). In the case of Cantor ternary set in Example 1.1, where c = 3, its 
dimension is dim.H{A) = dimB{A) = log 2/ log 3. □ 

Example 1.3. Let X C C is a triangle with vertices 0, a, and 1, where a is a complex 
number. Consider an iterated function system {X, F), where F = {/i, with fi{z) = az, 
and f2iz) = (1 — a)'z + a, where a satisfies \a\ < 1 and |1 — a| < 1, and z is the complex 
conjugate of z. (This case was proposed by De Rham [Dr].) 

If a = 1/2 -I- i\/3/6, then X is equilateral, and the attractor of this iterated function 
system is the Koch curve K, which can be constructed by replacing the middle 1/3 parts 
of each line segment of X with a new equilateral triangle, and doing this process with 
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iteration. In this case, |a| = |1 — a| = Ijsf^^ so their scahng numbers are Ai = A2 = l/\/3, 
and its similarity dimension is the unique solution of the equation 2(l/\/3)'' = 1- Thus 
dimij(iC) = dimB(ir) = log 4/ log 3. 

If \a — 1/2| = 1/2, e.g. a = 1/2 + i/2, then the attractor of {X,F) is the Peano curve 
P, which fully fill the space X. We have Ai = A2 = l/\/2, and solve the equation 
2(1/^2)*^ = 1, we obtain dimif(P) = dimB(P) = 2. □ 

2. Random Utility Functions 

We want to consider an agent in a large population who has a rational preference 
relation ^, but can be also affected by some psychological factors when making her choices. 
Concretely, we assume the decision state space is X, and the preference relation is a weak 
order ^ on X such that 

(i) either x ^ ?/ or y ^ x for all x, ?/ G X, 

(ii) x ^ y and y z implies x ^ z for all x,y,z E X. 

If 3; ^ y and not y ^ x, we say x)^ y. Then we have the following preference represen- 
tation result. 

Claim 2.1. If X is countable, or X is uncountably infinite but ^ is dense on X, which 
means that for any Z (Z X, when x,y & X\Z such that x y y, there is a z & Z satisfying 
X y z and z y y, then there exists a utility function u : X ^ R such that x ^ y if and 
only if u{x) > u{y) for all x,y E X. 

Let V{X) be the power set of X. A function C : V{X) — > P{X) is called a choice 

function, such that for any Z G 'P{X) we have C{Z) 7^ and C{Z) C Z. Now consider 
an arbitrary subset of X, say Y G V{X), if y G C{Y), then we should have u{y) > u{z) 
for all z gY, and if n(y) < u{z), then z ^Y. 

When we observe the agent w chooses y from a set of states we claim that u{y) > u{z) 
for any 2: G y is logically true. However, from the view of statistics, this inference is true 
only in a certain confidence level. So behind our observations, the utility of the agent w and 
also her preference is random for us, in the sense that the true utility values are adapted 
by noises. Such a problem was first considered by Manski [Ma], in which he introduced 
the method of maximum score to estimate the stochastic utility model of choice from the 
observation data. 

Define v{x) is the observed utility when the agent chooses x, then we have u{x) = 
v{x) + e(x), where e(x) is the introduced noise related to x. If v{x) is assumed to have a 
representation of linear form of m observable variables, say J(x) G M™, then v{x) = J{x)' j3., 
where P G is the parameter to be estimated. We have 

u{x) = J(x)'/3 + £(x), 

where the noises e(x) are i.i.d. across different alternatives x. The estimation $ is a 
function of the observed data J(x), where x belongs to a given sample, say 5 C X. 

Another way to deal with the difference between the utility associated with the observed 
choices and the true utilities is to think that the utility function of w is random for her, 
but the random utility function should have a fundamental unchanged part at any time. 
So w^s knowledge of her own preference ^ is probabilistic, and she uses a random utility 
function to make choices. The assumption that all the possible utility functions have a 
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common basic parts is set to capture the fact that she is not totally ambiguous with her 
own preference, but her knowledge is in a state between ignorance and certainty. 

Formally, suppose that there is a finite index set In with N > 2 elements, and the 
real- valued function Uj : X ^ R is a possible utility function for any i e In- Define a 
function k : X ^ M, we call it a fundamental utility function for the agent w, if there 
exists a contractive function fi'.X^X such that Ui = k o fi for all i £ In- 

So we have Ui{x) = fi(k{x)) for all i G In- Denote the image of the fundamental 
utility function k{x) on the domain X as k{X) C R. Then the random utility functions 
{ui : i G In} on X is equivalent with the function system {fi : i G In} on k{X). 

Now we introduce the time, and suppose the agent w makes her choices along the time 
domain Z_|_. At time t = 0, the fundamental utility function is ko{x) = k{x) for any x e X. 
At time t = 1, the fundamental utility function becomes ki{x) = fi{k{x)), which is also the 
utility function at time t = 0. In general, at time t = n, the fundamental utility function 
is kn{x) = fi{kn-i{x)) for any i & In and all n G N. This dynamics of fundamental utility 
functions can be thought of as a reasonable extension of the discounted utility function in 
the framework of a stationary time preference (^t)tez+ , where ^t=^ for all t G Z+ in our 
settings. 

This construction can help to generate an iterated function system {k{X),F), where 
F = {fi : i e In} with fi : k{X) — >■ k{X), and A; : X — >■ R is a given function. Since 
Ui{x) at time n is the fundamental utility function kn-\-i{x) at next time n + 1, the iterated 
function system {k{X),F) produces the values of both the fundamental utility and the 
actual utility. By the proposition 1.1, there exists an attractor A C k{X) such that 
A = Uie/^ fi{^)- Therefore, there exist some stable utilities under the iteration driven 
by F. 

If there is a probability measure n on In such that X^jg/j^ 7r(i) = 1, the iterated function 
system {k{X),F) becomes an iterated random function system {k{X),F,7r), in which we 
have the following random iterated process 

Un = fcn{Un-i), for all n G Z+, 

where Un is the random utility at time n, cj„ is the random index chosen from In at time 
n for all n G Z_|_, and U-i = k{x). Denote Kn as the random fundamental utility at time 
n, since i^n+i = Un, we have Kn+i = fa^{Kn) for all n G Z+ with Kq = k{x). 

Suppose fi{x) = piX, where < < 1 for all i E In, and pi ^ pj for any i ^ j, the 
iterated stochastic process is 

Un = ^nUn-1, for all n G Z+, 

where again Un is the random utility at time n, and f„ is an i.i.d. random variable with 
the support {pi : i G In} for all n G Z_|_. For any state x & X, at time n G Z_|_, the random 
utility is 

ran ra 

Un={ll Ci) k{x) = exp ( ^ log e^) k{x) = exp ( - ^ loglVeO) k{x), 

i=0 i=0 i=0 

where C/„ ^ almost surely, when n — >■ +oo. 

If fi{x) = px + Vi, where < p < 1 and rj > for all i G In, with 7^ rj for any i 7^ j, 
the iterated stochastic process is 

Un = pUn-l + En, for aU 71 G Z_|_, 
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where is again an i.i.d. random variable taking values from the set {rj 
n G Z+. At time n G Z_|., the random utility is 

n 

Un = p'%x) + ^p^-%, 
i=0 

in which p^k{x) vanishes when n — >■ +00, however the second part does not converge 
almost surely, simply as the new randomness e„ will enter into Un at each time n. Roughly 
speaking, in this iterated affine function system (k{X),F), the attractor will be a fractal 
for most possible values of the parameters p and n's. 

3. Optimal Growth under Uncertainty 

In this section, wc want to consider the nature of the steady state in a stochastic growth 
model. Consider an economy with a gross production function Y = F{K,L), where Y is 
the total production in the economy, K is the capital input, and L is the labor supply. 
Suppose this gross production function is a homogeneous function of degree 1, then we 
have Y/L = F{K/L, 1) = f{K/L). Define y = Y/L and k = K/L, we have the production 
function for a representative agent w in that economy y = f{k), which satisfies the classical 
hypothesis /(O) = 0, and f'{k) > 0, f"{k) < 0, for all k > 0, and the Inada conditions of 
limjfc_^o+ f'ik) = +00, and limfc^+00 f'{k) = 0. 

Suppose there is a stochastic factor ^ can affect the production, we write the new 
production function as y = f{k,^). In general, the shock ^ can be either multiplicative or 
additive with f{k), similar with the studies carried by Mitra, Montrucchio, and Privileggi 
[MMP], wc consider the case of multiplicative shocks. So we can write y = f{k, ^) = ^f{k), 
where ^ > is a random shock. Suppose ^ can take values in the set {Aj : i G In}, where 
In is a finite index set with N > 2, with the probability measure tt on In such that 

J2i&iN ■^(O ~ ^i''') > for all i G In- 

The sustainability of this economy is possible by introducing the consumption C and the 
investment I, both of whom have the same measurement with the production Y. Again 
we consider the consumption and the investment of the representative agent w, and hence 
we write her consumption as c = C/L, and her investment as i = I/L. 

After introducing the time domain Z_|_, we have the set of variables {yn,kn,imCn,^n} 
at each time n G Z+. In this economy without the government, we have the following 
system 

Vn ~ (hi in 
< in — , 
^ Vn = Cnfikn) 

where at time t = 0, ko = i_i is the initial capital, which produce yo = Cof{ko), and j/o 
will be divided into the consumption cq and the investment zq = fci. In general, at time 
t = n, kfi = in-i is the capital, and the production is yn = ^,nf{kn), which is the source of 
the consumption c„ and the next time capital kn^i = in, where ^„ is the random shock ^ 
for any n G so the process (^n)n is i.i.d. on the time domain Z_|_. 

Since in and kn+i are related by z„ = kn+i, and c„ and in are related hy Cn = yn — in 
given yn, we actually have the following three equivalent systems, which have the sets of 
variables {yn,kn+i,Cn}, {yn,Cn,Cn}, and {yn,in,Cn} respectively. 
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Suppose w has a stationary utility function n(c) with the properties that u'{c) > 0, 
u"{c) < 0, and limc^o+ ^'(c) = +oo, which means c„ > at each time n G Z4.. tw's time 
preference is captured by a discounting factor p G (0,1), so her utility along Z_|_ can be 
written as X^^g^^ p'^u{cn)-, if given her consumption flow (cq, ci, . . . , c„, . . . ). 

The steady states of this economic system are the solutions for the following dynamic 
optimization problem 

max Eq p^u{cn) 

S.t. kn+l = ^nfikji) — Cji 

fco > is given 

where is an i.i.d. random variable with a finite support {Aj : i G /at} and a probability 
measure vr, such that Xlie/jv ~ 

The optimal consumption flow of w is captured by the following Euler equation 

u'{Cn) = Et{p^n+lf'{kn+l)u'{Cn+l)). 

Noticing that kn+i = Dn — Cn-, we can rewrite the above equation as 

u'{Cn) = pfiVn - Cn)^t{^n+iu' (Cn+l)) . 

So for the optimal steady state along time, there exists a function : M ^ M such that 
c„ = ifiyn), and hence kn+i = Un — viVn)- Therefore, the optimal state of the production 
in this economy has the following dynamics 

Vn+l = Cn+lf{kn+l) = ^n+lf [Vn - ViVn)) = in+li^iVn), 

where ip{yn) = /{Vn — v{yn)), and ^„ is chosen randomly from {A^ : i G In} according to 
the probability measure tt on 1^. 

Suppose there exists an invariant support set (related with the 7r-balanced measure for 
this Markov chain) for the stochastic process yn+i = Cn+iV'(yn)j say an compact interval 
Xy located within M_|_. Define gi{y) ■= Xii/j{y) for all i G In- Then we have an iterated 
random function system (Xy,G, tt), where G = {gi : i E In} with gi : Xy Xy for all 
i & In, and tt is a probability measure on In- 

We also have the dynamics of the optimal states of the capital in this economy 

kn+l =yn- ViVn) = ^nf{kn) - ^{Cnf{kn)) = 4>{kn,in), 

where (l){kn,^n) = ^nfikn) - f{^nf{kn))- Define hi{k) := (t){k,Xi) for all i G In- Let the 
invariant support interval for fc^+i = (l>{knj Cn) be Xk x {Aj : i G In}- We have an iterated 
function system (Xk, H), where H = {hi : i £ In} with hi : Xk X^ for all i £ In- 

Consider a simple case, where f{k) = and w(c) — log c. Suppose the shock is 
an i.i.d. Bernoulli process, such that ^„ = Aa with probability 7r(a) = q and = A;, with 
probability 7r{b) = 1 — q, here In = {a,,h}. Wc assume l/Af, > Aa > 1 > Ab > 0, so the 
index a stands for a positive shock, and the b stands for a negative shock, but the negative 
shock will not make the economy vanish, and at the same time the positive shock will not 
make it too expansive. 

Then at the steady states of its optimal growth, we have Cn = ^{yn) = (1 — p/'^)yn, so 
kn+i = yn-Cn = pyn/^- We know yn = Cnf{kn) = ^n\/K, SO kn+l = Define 
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Kn = log kn, then we have 

Kn+l = Kn/S + (log^„ + log /9 - logS). 

So when ^„ takes the value Xa or Ab, we have two affine functions, namely lai^) = k/3 + 
(log Xa + log p - log 3) and lb{K) = k/3 + (log A;, + log /> - log 3). 

Suppose the invariant support interval of la{n) and Ibin) is [a, then we have (3/3 + 
(log Xa + log p — log 3) = /3, and q;/3 + (log Af, + log p — log 3) = a. After some calculations, 
we obtain a = 3(log A5 + logp — log3)/2, and /3 = 3(logAa + logp — log3)/2. So there 
exists an iterated function system ([a,f3],L), where L = {la,lb} with li : [a, P] — >■ [a,P] for 
i = a,b. 

The difference between these two /(/«)-intercepts are logAa — logA^ = log(Aa/Ab), and 

(3 — a = 3(log Xa — log Aj) /2 = 3 log(Aa/A5) /2. So after some transformations, the function 
system lain) and lb{n) is equivalent with the pair of functions QaiT) = r/3 and ^(t) = 
r/3 + 2/3, where t = k/{/3 — a) + d, where d is a constant. So the iterated function system 
{[a,/3],L) is equivalent with the iterated function system ([0,1], Q), where Q = {qa,Qb} 
with Qi : [0, 1] —7- [0, 1] for i = a,b. From Example 1.1, we know the attractor of ([0, 1], Q) 
is the Cantor ternary set, so the attractor of ([a, (3],L) is also a Cantor set, and hence the 
dynamics of the steady states of Kt and also kt = exp(«;t) is in a chaotic situation. 

4. Summary 

In this work, we study the random utility function and the optimal stochastic growth 
based on the formal developed theory of iterated function systems. Since the real economic 
world is very complex, we have to refine again and again our model settings and theoretical 
hypothesis to have a considerably thinkable framework. However, I wish these frameworks 
were able to capture some basic natures of the economic phenomena. 
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